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Abstract. About half of all known stellar systems with Sun-like stars consist 
of two or more stars, significantly affecting the orbital stability of any planet in 
these systems. Here we study the onset of instability for an Earth- type planet 
that is part of a binary system. Our investigation makes use of previous analyt- 
ical work allowing to describe the permissible region of planetary motion. This 
allows us to establish a criterion for the orbital stability of planets that may be 
useful in the context of future observational and theoretical studies. 



Introduction 



Observational evidence for the e xistence of planets iii stellar binary (and higher 
order) systems has be e n give n bv lPatience et all (120021') .lEggenbereer et alJ (|2004l ') , 
Eggenberger Udrvl ( 20071 ). and others. Eggenberger Udrv presented data 



for more than thirty systems, mostly wide binaries, as well as several triple star 
systems, with separation distances as close as 20 AU (GJ 86). These observations 
are consistent with the finding that binary (and hig her order) systems occur in 
high freque i icy in the local Galact i c neighborhood (iDuquennov &: Mayor 1991 



iLadal l2006l : iRaghavan eFaD l2006l : iBonavita Desideral l2007l ) . The fact that 
planets in binary systems are now considered to be relatively common is also 
implied by the recent detection of debris disks in vario us main-sequence s tellar 
binary systems using the Spitzer Space Telescope (e.g.. iTrilling et 311120071 ). 

In the last few decades, significant progress has been made in the study 
of stability of planetary orbits in stellar binary systems. Most of these studies 
focused on S-type systems, where the planet is orbiting o ne of the sta r s with 
the second star to be considered a perturbator. Recently, iDavid et al.l ( 20031 ) 
investigated the orbital stability of an Earth-mass planet around a solar-mass 
star in the presence of a companion star and determined the planet's ejection 
time for systems with a vari ety of orbit a l eccentricities and se mimajor axes. 

In our previous work (jStuitl Il995l : iMusielak et"all l2005l ). we studied the 
stability of both S-type and P-type orbits in stellar binary systems, and deduced 
orbital stability limits for planets. These limits were found to depend on the 
mass ratio b e tween the stellar compo n ents. This topic was recently revisited by 
Cuntz et al.l ( 20071 ) and lEberle et al. ( 2007 ) , who used th e concept of Jacobi's 



integral and Jacobi's constant ( SzebehehT 19671 : Rov 20051 ) to deduce stringent 
criteria for the stability of planetary orbits in binary systems for the special 
case of the "coplanar circular restricted three-body problem" . In this paper, we 
present case studies of planetary orbital stability for different stellar mass ratios 
and different initial planetary distances from its host star. 
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Figure 1. Model simulations for /x = 0.2 and 0.5, and different values of po- 

Each panel shows the primary (large dot) and secondary {small dot) star, the 
planetary orbit (solid line), the five Lagrange points, and the "zero velocity 
contours" (dash-dotted lines). The Lagrange points are denoted as L2, LI, L3, 
respectively, from left to right along the line connecting the two stars, and as 
L4 (top) and L5 (bottom) apart from this line. For /x = 0.2, the critical values 

p'^\ p^q \ and Po^^ are given as 0.353, 0.420, and 0.692, respectively, and for 
p, = 0.5, they are given as 0.251, 0.442, and 0.442, respectively. 
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2. Methods and results 

In the so-called coplanar circular restricted three-body problem the two stars are 
assumed to orbit each other in circles and their masses are much larger than that 
of the planet. In our case, it is assumed that the planetary mass is 1 x 10~^ of 
the mass of the star it orbits; also note that the planetary motion is constrained 
to the orbital plane of the two stars. In addition, it is assumed that the initial 
velocity of the planet is set for an initially circular orbit, and that it is in the 
same direction as the orbital velocity of its host star. This star shall be the more 
massive of the two stars. Furthermore, the starting position of the planet is to 
the right of its host star along the line joining the binary components (3 o'clock 
position). The mass ratio /x of the two stars is defined as // = M2/M with 
M = Ml + M2 , where Mi and M2 are the masses of the primary and secondary 
star, respectively. Additionally, po denotes the planet's relative initial distance 
po = Rq/D from its host star, with D as distance between the two stars and Rq 
the initial planetary distance from the primary star. 

In the following, we illustrate the transition from stability to instability by 
progressively increasing the value of po for binary systems with a fixed mass 
ratio /U, given as p, = 0.2 and 0.5, respectively. For both p = 0.2 and 0.5, we 
present the resulting planetary orbits for four different values of po (see Fig. 1). 
Each panel shows the primary (larger dot) and secondary (smaller dot) star, the 
planetary orbit (solid line), as well as the zero velocity contours (dash-dotted 
lines). The upper four panels of Fig. 1 refer to the stellar mass ratio p = 0.2, 
whereas the lower four panels refer to p = 0.5, 

Let us first focus on the case studies for p = 0.2 and po = 0.25 and 0.323. 

Both values of po are smaller than the critical value of p'q^ = 0.353, indicating 
that the planetary orbits are stable. The fact that we restricted the time of 
the simulation to 50 yrs is inconsequential owing to the fact that the orbital 
stability of the planet is guaranteed by the analytical properties of the system, 

namely po < pg^^; see ICuntz et al.l ( 20071 ) and lEberle et al.l ( 20071 ) for a more 



extended discussion. Also note that the zero velocity contour changes between 
the two panels due to the increase in pQ by getting closer to the Lagrange point 
LI, although its topology remains unaltered. Moreover, we show two cases of 
unstable orbits by choosing po = 0.41 and 0.51, respectively. Since both values 

exceed p'^\ the zero velocity contour opens at LI, providing the possibility 

for the planet to be captured by the secondary star. For pQ = 0.51, the zero 

(2) (3) 
velocity contour even opens at L2 because of po > Po • ^^^^ Po > Po (i^ot 

shown here), the contour would even open at L3, providing a further type of 
opportunity for the planet to escape from the binary system. 

The results for po = 0.41 show that the planetary orbit is unstable but still 
remains within the sphere of gravitational influence of the primary star. The 
situation is different for pQ = 0.51 where the planet reaches the secondary star 
only after a few irregular orbits about the primary star have been completed. 
A more detailed analysis shows that the planet first encountered the secondary 
star after 39.6 yrs at a minimal distance of at most 0.1 AU. A second even closer 
encounter occurred after 41.2 yrs, when the simulation was stopped because the 
planet entered the Roche limit of the secondary star. The general behavior of 
the model with po = 0.51 is due to the fact that po exceeds both pg^^ and pg^\ 
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which in principle allows the planet to escape from the binary system through 
the L2 point. 

Our results for /i = 0.2 demonstrate different cases of orbital stability and 
instability. Similar results are obtained for /i = 0.5, albeit quantitative differ- 
ences due to the different values of //, Pq^\ Pq \ and Pq\ We find again that 

orbital stability is obtained if pQ < Pq^^ , whereas for larger values of po instability 
is expected to emerge. Highly unstable cases are found for po = 0.3 and 0.45. 



3. Conclusions 

For the special case of the "coplanar circular restricted three-body problem", 
we applied stringent mathematical criteria that allow to precisely determine 
whether a planetary orbit in a stellar binary system is stable or unstable. This 
is accomplished by comparing the planet's relative initial distance po to the 

critical values Pf^\ Pq'\ and p^q \ defined for a fixed stellar mass ratio p. An 
adequate way of demonstrating this different type of behavior is the assessment 
of the topology of the zero velocity contour, given by the p and /?o values of 
the system. In this case, planetary orbital stability is obtained if the contour 
is completely closed around the primary star and planet. All numerical case 
studies show a behavior consistent with this theoretical prediction. Note how- 
ever that for generalized binary systems, other meth ods are required to deter- 
mine the long-ter m stability of planetary orbits (e.g., Holman Wiegert 19991 : 



David et al. 20031 ) . Important applications of our work include contesting nu- 



merically deduced stability limits for cases where analytically deduced results 
exist. 
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